ON 

o 
o 

> 

o 



Two-flavor lattice QCD study of the axial charges of N(1535) and N(1650) 

Toru T. Takahashi a and Teiji Kunihiro b 

a Yukawa Institute for Theoretical Physics, Kyoto University, 
Kitashirakawa-Oiwakecho, Kyoto 606-8502, Japan 

b Department of Physics, Kyoto University, 
Kitashirakawa-Oiwakecho, Kyoto 606-8502, Japan 



Oh! 
43 



> 

m 
(N 



On 
O 



We show the first lattice QCD results on the axial charges of iV*(1535) and iV*(1650). The measurements 

are performed with two flavors of dynamical quarks employing the renormalization-group improved gauge action 
at /3=1.95 and the mean-field improved clover quark action with the hopping parameters, k=0.1375, 0.1390 
and 0.1400. We construct 2x2 correlation matrices and diagonalize them so that the signals of iV*(1535) and 
iV*(1650) are properly separated. Wraparound contributions in the correlator, which can be another source of 
signal contaminations, are eliminated by imposing the Dirichlet boundary condition in the temporal direction. 
The axial charge of iV*(1535) is found to take small values as ~ O(0. 1), whereas that of iV*(1650) is 

approximately 0.5, which is almost independent of quark masses and consistent with the predictions by the naive 
nonrelativistic quark model. 



1. Introduction 

Crural symmetry together with its spontaneous 
breaking is one of the key ingredients in the 
low-energy hadron or nuclear physics. Due to 
its spontaneous breaking, up and down quarks, 
whose current masses are of the order of a few 
MeV, acquire the large constituent masses of a 
few hundreds MeV, and are consequently respon- 
sible for about 99% of mass of the nucleon and 
hence that of our world. One would thus consider 
that chiral condensate (ipip), the order parame- 
ter of the chiral phase transition, plays an essen- 
tial role in the hadron-mass genesis in the light 
quark sector. On the other hand, chiral symme- 
try gets restored in systems where hard external 
energy scales such as high-momentum transfer, 
temperature(T), baryon density and so on exist, 
owing to the asymptotic freedom of QCD. Then, 
several questions may arise: Are all hadronic 
modes massless in such systems? Can hadrons 
be massive even without non- vanishing chiral con- 
densate? 

An interesting possibility was suggested some 
years ago by DeTar and Kunihiro [1] , who showed 
that nucleons can be massive even without the 



help of chiral condensate due to the possible chi- 
rally invariant mass terms, which give degener- 
ated finite masses to the members in the chiral 
multiplct (a nucleon and its parity partner) even 
when chiral condensate is set to zero: In order 
to demonstrate this possibility for a fmite-T case, 
they introduced a linear sigma model which has 
a nontrivial chiral structure in the baryon sector 
and a mass-generation mechanism essentially dif- 
ferent from that by the spontaneous chiral sym- 
metry breaking. Interestingly enough, their chi- 
ral doublet model has recently become a source 
of debate as a possible scenario of observed parity 
doubling in excited baryons 2 3 4 5 6 7J, although 
their original work [1] was intended for an appli- 
cation to fmite-T systems. 

It is thus an intriguing problem to clarify the 
chiral structure of excited baryons in the light 
quark sector beyond model considerations. One 
of the key observables which are sensitive to 
the chiral structure of the baryon sector is ax- 
ial charges [I]. The axial charge of a nucleon N 
is encoded in the three-point function 

(N\A^\N) = u^-[ W g A (q 2 ) + q^ 5 h A (q 2 )]u.(l) 

Here, A a = Qj^Js^ttQ is the isovector axial cur- 
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rent. The axial charge qa is defined by gA (q 2 ) 
with the vanishing transferred momentum q 2 = 0. 
It is a well-known fact that the axial charge N 
of (940) is 1.26. Though the axial charges in the 
chiral broken phase can be freely adjusted with 
higher-dimensional possible terms and cannot be 
the crucial clues for the chiral structure [314] , they 
would reflect the internal structure of baryons and 
would play an important role in the clarification 
of the low-energy hadron dynamics. 

In this paper, we show the first unquenched 
lattice QCD study [8] of the axial charges g 1 ^ N 
of 7V*(1535) and 7V*(1650), two lowest nucleon 
resonances in the negative parity channel. We 
employ 16 3 x 32 lattice with two flavors of dy- 
namical quarks, generated by CP-PACS collab- 
oration [9] with the renormalization-group im- 
proved gauge action and the mean-field improved 
clover quark action. We choose the gauge config- 
urations at (3 — 1.95 with the clover coefficient 
csw = 1-530, whose lattice spacing a is deter- 
mined as 0.1555(17) fm. We perform measure- 
ments with 590, 680, and 680 gauge configura- 
tions with three different hopping parameters for 
sea and valence quarks, k sc& , K va i = 0.1375, 0.1390 
and 0.1400, which correspond to quark masses of 
~ 150, 100, 65 MeV and the related ir-p mass 
ratios are m PS /m v = 0.804(1), 0.752(1) and 
0.690(1), respectively. Statistical errors are esti- 
mated by the jackknife method with the bin size 
of 10 configurations. 

2. Formulation 

Our main concern is the axial charges of the 
negative-parity nucleon resonances AT* (1535) and 
iV*(1650) in i channel. We then have to con- 
struct an optimal operator which dominantly cou- 
ples to TV* (1535) or TV* (1650). We employ the 
following two independent nucleon fields, 

N x {x) = e ahc u a (x)(u b (x)C 75 d c (x)) (2) 

and 

N 2 (x) = e ahcl5 u a (x)(u b (x)Cd c (x)), (3) 

in order to construct 2x2 correlation matri- 
ces and to separate signals of AT* (1535) and 
AT*(1650). Here, u(x) and d(x) are Dirac spinors 



for u- and d- quark, respectively, and a, 6, c denote 
the color indices. Even after the successful signal 
separations, there still remain several signal con- 
taminations: Signal contaminations by scattering 
states and wraparound effects. 

Due to the unquenched gauge configurations, 
the negative parity nucleon states could decay 
to 7r and N, and their scattering states could 
inevitably get into the spectrum. The sum of 
the pion mass M„ and the nucleon mass Mn is 
however in our setups heavier than the masses 
of the lowest two states (would-be AT* (1535) 
and AT* (1650)) in the negative parity channel. 
We then do not suffer from any scattering-state 
signals. The other possible contamination is 
wraparound effects [lOj . Since we perform un- 
quenched calculations, the excited nucleon N* 
can decay into N and n, and even when we have 
no scattering state |iV+7r), we could have another 
type of "scattering states" . The baryonic correla- 
tor (N*(t sn k)N*(t SIC )) can still accommodate, for 
example, the following term. 

(7r\N*(t snk )\N)(N\N*(t SIC )\n) 

X g-- E N(t S nk-tsrc) x g--E^(JV t -t 31lk + t src ) /^\ 

Here, N t denotes the temporal extent of a lattice. 
Such a term is quite problematic and mimic a 
fake plateau at En — E n in the effective mass plot 
because it behaves as ~ e~( EN ~ E ^ l - ta,ik ~ tsrc \ In 
order to eliminate such contributions, we impose 
the Dirichlet condition on the temporal boundary 
for valence quarks, which prevents valence quarks 
from going over the boundary. (Wraparound 
effects can be found even in quenched calcula- 
tions [TO]/) 

We employ zero-momentum-projected point- 
type operators for the sinks, and employ wall- 
type operators in the Coulomb gauge for the 
sources. After we diagonalize the 2x2 corre- 
lation matrices, we can construct optimized op- 
erators for N(1535) and N(1650) states. Once 
we construct optimized operators, we can easily 
compute the (non-renormalized) vector and ax- 
ial charges ^ or ^ ne positive- and negative- 
parity nucleons via three-point functions with the 
so-called sequential-source method [TT]. In prac- 
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tice, we evaluate g^ lat \t) denned as 



Tr T A (B(t 
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Tr r^(S(i 8nk )B(i src )) 



(5) 



and extract g^' lat ' by the fit g^' lat ' 



±[lat]/,N . 
9 A (*) m 



the plateau region. denotes the (optimized) 
interpolating field for nucleons. Ta is 7^75 "^ 74 , 
and Jfi(t) is an axial vector current inserted at t. 

We show in Fig. [T] the non-renormalized axial 
charge g^°' lat ' (t) for N* (1535) as a function of the 
current insertion time t. They are rather stable 
around t STC < t < t sn k. 
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Figure 1. The non-renormalized axial charge of 
TV* (1535), 3^°' lat '(i), as a function of the current 
insertion time t. 

We finally reach the renormalizcd charges 
g\ = ZaQ^^ with the prefactors Za = 
2ku Za (l + bA^^j , which are estimated with 
the values listed in Ref. [9]. 

3. Results 

We first take a stock of the axial charge g^~ of 
the ground-state positive-parity nucleon, which 
is well known and can be the references. The 
axial charge of the positive parity nucleon 
is shown in Fig. [5] as open squares. The axial 
charge of the positive parity nucleon shows little 



Figure 2. The renormalized axial charges of the 
positive- and the negative-parity nucleons are 
plotted as a function of the squared pion mass 
m^.. The solid line is drawn at gA = 1.26 and the 
dashed line is drawn at gA = 0. 



quark-mass dependence, and they lie around the 
experimental value 1.26. 

We finally show the axial charges of the 
negative-parity nucleon resonances in Fig. [5] One 
finds at a glance that the axial charge <n - of 
AT* (1535) takes quite small value, as g^ ~ 
O(0. 1) and that even the sign is quark-mass de- 
pendent. While the wavy behavior might come 
from the sensitiveness of gT to quark masses, 
this behavior may indicate that g ^ is rather 
consistent with zero. These small values are 
not the consequence of the cancellation between 
u- and d-quark contributions. The u- and d- 
quark contributions to g^ are in fact individ- 
ually small [5]. On the other hand, the axial 
charge of iV*(1650) is found to be about 
0.55, which has almost no quark-mass depen- 
dence. The striking feature is that these axial 
charges, g^~ ~ and g 1 ^ ~ 0.55, are consistent 
with naive nonrelativistic quark model calcula- 
tions [12113] . g\ = — I and g^ - = |. Such values 
are obtained if we assume that the wave functions 
of 7V*(1535) and iV*(1650) are \l=l,S=%) and 
\l — 1,S — |) neglecting the possible state mix- 
ing. (Here, I denotes the orbital angular momen- 
tum and S the total spin.) 
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In the chiral doublet model [112] . the small 
g 1 ^ N is realized when the system is decoupled 
from the chiral condensate (i/jip). The small 
g^~ of iV*(1535) then does not contradict with 
the possible and attempting scenario, the chi- 
ral restoration scenario in excited hadrons [2]. If 
this scenario is the case, the origin of mass of 
iV*(1535) (or excited nucleons) is essentially dif- 
ferent from that of the positive-parity ground- 
state nucleon 7V(940), which mainly arises from 
the spontaneous chiral symmetry breaking. How- 
ever, the non-vanishing axial charge of TV* (1650) 
unfortunately gives rise to doubts about the sce- 
nario. 

4. Conclusions 

We have performed the first lattice QCD study 
of the axial charges g 1 ^ N of iV*(1535) and 
iV*(1650), with two flavors of dynamical quarks 
employing the renormalization-group improved 
gauge action at /3=1.95 and the mean- field im- 
proved clover quark action with the hopping pa- 
rameters, k=0.1375, 0.1390 and 0.1400. We have 
found the small axial charge g°^~ of TV* (1535), 
whose absolute value seems less than 0.2 and 
which is almost independent of quark mass, 
whereas the axial charge g\~ of TV* (1650) is found 
to be about 0.55. These values are consistent 
with naive nonrelativistic quark model predic- 
tions, and could not be the favorable evidences for 
the chiral restoration scenario in (low-lying) ex- 
cited hadrons. Further investigations on the axial 
charges of TV* (1535) or other excited baryons will 
cast light on the chiral structure of the low-energy 
hadron dynamics and on where hadronic masses 
come from. 
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